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We present a method for calculating transfer matrices for the g-state Potts 
model partition functions Z{G,q,v), for arbitrary q and temperature variable 
V, on cyclic and Mobius strip graphs G of the square (sq), triangular (tri), 
and honeycomb (he) lattices of width Ly vertices and of arbitrarily great 
length Lx vertices. For the cyclic case we express the partition function as 
Z{A,Ly X Lx,q,v) = EjloC^'^^^^[(^^,A,L„d)'"], where A denotes lattice type, 
c^'^^ are specified polynomials of degree d in q, Tz^A,Ly,d is the transfer matrix in 
the degree-d subspace, and m = {Lx/2) for A = sq, tri (he), respectively. 
An analogous formula is given for Mobius strips. We exhibit a method for 
calculating Tz^\.Ly.d for arbitrary Ly. Explicit results for arbitrary Ly are 
given for Tz^A,Ly,d with d = Ly and d = Ly — 1. In particular, we find 
very simple formulas the determinant det{Tz^A,Ly,d), and trace Tr(Tz^A,Ly)- 
Corresponding results are given for the equivalent Tutte polynomials for these 
lattice strips and illustrative examples are included. We also present formulas 
for self-dual cyclic strips of the square lattice. 
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I. INTRODUCTION 



The g-state Potts model has served as a valuable model for the study of phase transitions 
and critical phenomena [1,2]. On a lattice, or, more generally, on a (connected) graph G, at 
temperature T, this model is defined by the partition function 

Z(G,g,^) = 5]e-^« (1.1) 

{CT„} 

with the (zero-field) Hamiltonian 

n = -jJ2Sa.a, (1.2) 

{ij) 

where = 1, ...,q arc the spin variables on each vertex (site) i E G\ (3 = (kBT)^^; and (ij) 
denotes pairs of adjacent vertices. The graph G = G{V, E) is defined by its vertex set V and 
its edge set E] we denote the number of vertices of C as n = n{G) = \V\ and the number of 
edges of G as e{G) = \E\. We use the notation 

K ^(3J , v^e^ -1 (1.3) 

so that the physical ranges are v > ior the Potts ferromagnet, and — 1 < ^; < for Potts 
antiferromagnet, corresponding to < T < oo. One defines the (reduced) free energy per 
site / = — where F is the actual free energy, via /({G}, q, v) — hm„_»oo ln[Z(G, q, vy^'"], 
where we use the symbol {G} to denote the formal limit lim„_^oo G for a given family of 
graphs. In the present context, this n — > oo limit corresponds to the limit of infinite length 
for a strip graph of the square lattice of fixed width and some prescribed boundary conditions. 

In this paper we shall present transfer matrices for the g-state Potts model partition 
functions Z{G,q,v), for arbitrary q and temperature variable v, on cyclic and Mobius strip 
graphs G of the square, triangular, and honeycomb lattices of width Ly vertices and of 
arbitrarily great length vertices. We label the lattice type as A and abbreviate the three 
respective types as sq, tri, and he. Each strip involves a longitudinal repetition of m copies 
of a particular subgraph. For the square-lattice strips, this is a column of squares. It is 
convenient to represent the strip of the triangular lattice as obtained from the corresponding 
strip of the square lattice with additional diagonal edges connecting, say, the upper-left to 
lower-right vertices in each square. In both these cases, the length is = m vertices. We 
represent the strip of the honeycomb lattice in the form of bricks oriented horizontally. In 
this case, since there are two vertices in 1-1 correspondence with each horizontal side of a 
brick, — 2m vertices. Summarizing for all of three lattices, the relation between the 
number of vertices and the number of repeated copies is 
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T _ ( m if A = sq OT tri or Go /i a\ 

~ I 2m if A = /ic ^ ' ' 

Here Gd is the cyclic self-dual strip of the square lattice, to be discussed further below. For 
the cyclic case the partition function has the general form [11,12] 

Ly nz(A,Ly,d) 

Z{A, Ly X L„ eye, q,v)^J2 ^^'^ E i^WAjr (1-5) 

d=0 j=l 

In terms of transfer matrices, this can be written as 

Ly 

Z(A, Ly X L„ eye, q,v)^Yl c^'^Tr[{Tz,A,Ly,dr] (1-6) 

d=0 

where c^^^ are polynomials of degree d in g defined below, and the transfer matrices in the 
degree-d subspace Tz,K,Ly,d and their eigenvalues Xz,A,Ly,d,j, j = 1, ...,nz(A, Ly,d) are inde- 
pendent of the length of the strip length L^.. Here we shall construct an analogous formula 
for Mobius strips. We exhibit our method for calculating T^.A.Ly.d for arbitrary Ly. Ex- 
plicit results for arbitrary Ly are given for (i) Tz^K^Ly,d with d = Ly and d — Ly — 1, (ii) 
the determinant det(Tz^A^Ly,d), and (iii) the trace Tr{Tz^A,Ly)- Corresponding results are 
given for the equivalent Tutte polynomials for these lattice strips and illustrative exam- 
ples are included. We have calculated the transfer matrices up to widths Ly = 5 for the 
square, triangular, and honeycomb lattices and Ly = 4 for the cyclic self-dual strip of the 
square lattice. Since the dimensions of these matrices increase rapidly with strip width (e.g., 
dim(Tz^A,Ly,d) = 14, 28, 20, 7 for Ly = 4 and < < 3), it is not feasible to present many of 
the explicit results here; instead, we concentrate on general methods and results that hold 
for arbitrary Ly. 

Various special cases of the Potts model partition function Z{G, q, v) are of interest. For 
example, if one considers the case of antiferromagnetic spin-spin coupling, J < and takes 
the temperature to zero, so that K — —oo and v — —1, then 

Z(G,q,-l)^P(G,q) (1.7) 

where P{G, q) is the chromatic polynomial (in q) expressing the number of ways of coloring 
the vertices of the graph G with q colors such that no two adjacent vertices have the same 
color [8,9]. 

We recall some previous related work. The partition function Z{G,q,v) for the Potts 
model was calculated for arbitrary q and v on strips of the lattice A with cyclic or Mobius 
longitudinal boundary conditions was calculated for (i) A — sq, Ly — 2 in [11] and the 
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nature of the coefficients c^*^^ and subspace dimensions n^(A, Ly, d) given in [12], (ii) A = sg, 
= 3 in [13], (iii) A = tri, Lj^ = 3 in [14], (iv) for A = /ic in [15], (v) for cyclic self-dual 
strips of the square lattice with Lj^ = 1, 2, 3 in [16-18], and (vi) for = 2, 3 on the square- 
lattice with next-nearest-neighbor spin-spin couplings in [19,20]. Some general structural 
properties for lattice strips with arbitrary width were given in [12,17,21]. Matrix methods 
for calculating chromatic polynomials were developed and used in [22,23], [24,25] and more 
recently in [26]- [28]. Ref. [29] developed transfer matrix methods for both Z{G,q,v) and 
the special case v = —1 oi chromatic polynomials on strips of the square lattice with free 
longitudinal boundary conditions and used them to calculate the latter polynomials for a 
large variety of widths. These have been termed transfer matrices in the Fortuin-Kasteleyn 
representation (see eq. (1.8) below). These methods were apphed to calculate the full Potts 
model partition function for strips of the square and triangular lattices with free boundary 
conditions and a number of widths in Refs. [30,31]. A number of calculations have been 
done for the special case of chromatic polynomials of lattice strip graphs; we do not review 
these here but refer the reader to the references in [34] . In a different direction, we mention 
calculations for arbitrary q and v on rectangular lattice patches for the purpose of studying 
distributions of zeros [32,33]; we also do not review these studies here since our present focus 
is calculations for strips of arbitrarily great length. 

There are several motivations for presenting these transfer matrix methods for calculat- 
ing Potts model partition functions on cyclic and Mobius lattice strips and general results 
that we have obtained for lattice strips of arbitrary width as well as length. Clearly, new 
exact results on the Potts model are of value in their own right. The transfer matrices are 
also a convenient way to calculate partition functions. In explaining this, one should note 
that a very compact way of expressing the relevant information is in the form of generat- 
ing functions for Z{G,q,v) or in the form where the eigenvalues are given as roots of the 
characteristic polynomials. For a particular degree-d subspace, it requires nz{A, Ly,d) co- 
efficients to specify these characteristic polynomials. In contrast, it requires nz{A, Ly,dy 
entries (some of which may be the same) to specify the transfer matrix Tz^\^Ly,d- However, 
although the transfer matrix is not the most compact way of presenting the necessary in- 
formation to specify Z{G,q,v), it is, nevertheless, as convenient as the generating function 
method for calculating the partition function since one only needs to compute traces of pow- 
ers of the matrices. Moreover, except for the lowest few values of the strip width Ly, the 
characteristic polynomials are often of fifth order or higher, so that it is not possible to ob- 
tain explicit algebraic solutions for the eigenvalues. In our previous work we have made use 
of the theorem on symmetric polynomial functions of roots of an algebraic equation [36,37], 
which states that such functions can be expressed (via Newton identities) in terms of the 
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coefficients of tlie algebraic equation. Tlie particular type of symmetric polynomial function 
of the roots that is relevant for the Potts model partition functions on cyclic strips is the sum 
of m'th powers of these roots. Although the theorem on symmetric polynomial functions 
of roots of algebraic equations guarantees that these sums of m'th powers are expressible 
in terms of the coefficients of the equations, which arc polynomials in q and f , it does not 
imply that they are particularly simple, and, indeed, the Potts model partition functions 
of moderately long lattice strips arc quite lengthy expressions. Since the determinant and 
trace of a matrix are, respectively, the product and sum of the eigenvalues, they are both 
symmetric polynomial function of these eigenvalues (roots of the characteristic polynomial). 
The well-known expressions for the determinant and trace of a matrix in terms of coeffi- 
cients of the characteristic polynomial are examples of the above-mentioned theorem. The 
very simple formulas that we have obtained for determinants and traces of transfer matrices 
provide further motivation for the present exposition. 

Let G' — (y, E') be a spanning subgraph of G, i.e. a subgraph having the same vertex 
set V and an edge set E' C E. Z{G, q, v) can be written as the sum [3,4] 

Z{G,q,v)= /•(^'^^^l'"'' (1.8) 

G'CG 

where k{G') denotes the number of connected components of G' . Since we only consider 
connected graphs G, we have k{G) = 1. The formula (1.8) enables one to generalize q from 
Z+ to ]R+ for physical ferromagnetic v. More generally, cq. (1.8) allows one to generalize 
both q and v to complex values, as is necessary when studying zeros of the partition function 
in the complex q and v planes. 

The Potts model partition function Z{G, q, v) is equivalent to an object of considerable 
current interest in mathematical graph theory, the Tutte polynomial, T{G,x,y), given by 
[5]- [7] 

TiG,x,y) =Y{x- 1)^(«')-MG)(^ _ i)c(G') (19) 

G'CG 

where c(G') = \E'\ + k{G') — \V\ is the number of independent circuits in G' . Now let 

x = l + y^v + l (1.10) 

V 

so that 

q^{x-l){y-l) . (1.11) 

Then the equivalence between the Potts model partition function and the Tutte polynomial 
for a graph G is 
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Z{G, q, v)^{x- - X, y) . (1.12) 

Given this equivalence, we can express results either in Potts or Tutte form. We will use 
both, since each has its own particular advantages. The Potts model form, involving the 
variables q and v is convenient for physical applications, since q specifies the number of states 
and determines the universality class of the transition, and v is the temperature variable. 
The Tutte form has the advantage that many expressions are simpler when written in terms 
of the Tutte variables x and y. 

We recall that for a planar graph G — {V,E), one defines the (planar) dual graph G* 
as the graph obtained by replacing each vertex (face) of G by a face (vertex) of G* and 
connecting the vertices of the resultant G* by edges. The graph G is self-dual if and only 
ii G — G*. For a planar graph Gpi, it is evident from the definition (1.9) that the Tutte 
polynomial satisfies 

T{Gpi,x,y)=T{G;i,y,x) . (1.13) 

Equivalently, 

Z{G,i, q, v) = ^<«-')g-^(«^')Z(G;„ q, ^) . (1.14) 
The coefficients in eq. (1.5) are 

c^'^ = u,M"V'^) = Y.{-iy{^'^~\'-' (1-15) 

with Un{x) being the Chebyshev polynomial of the second kind. The first few of these 
coefficients are c*^°^ = 1, S-^^ = q — l, c*^^^ = — 3g + 1, and c^^^ = q^ — bq^ + 6g — 1. The c^'^^'s 
play a role analogous to multiplicities of eigenvalues \z,A.Ly,d,j, although this identification is 
formal, since c^'^'^ may be zero or negative for the physical values q = 1, 2, 3 [12]. For 5 > 4, 
c*^*^^ is positive-definite. From (1.12), one can write the Tutte polynomial as 

Ly nz{A,Ly,d) 

r(A, Ly X L„ cyc.,x, y) = —— c^'^ E (^^A^^dj)" (1-16) 

d=0 j=l 

where m is given in terms of by eq. (1.4) and it is convenient to factor out a factor of 
l/{x — 1). (This factor is always cancelled, since the Tutte polynomial is a polynomial in x 
as well as y.) In terms of transfer matrices, 

T{A, Ly X L„ eye, x, y) = —— J] cWTr[(rT,A,L„d)™] • (1.17) 

d=Q 
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The equivalence of eq. (1-17) and (1.6) to (1.16) and (1.5), respectively, relies upon the 
theorem that an arbitrary square matrix T can be put into (upper or lower) triangular form 
(tf) by a (unitary) similarity transformation [35] 

Ijrpjj-i ^ rp^^ _ (118) 

For dcfiniteness, we consider the upper triangular form, for which {Ttf)ij = if i > j. 
A matrix in triangular form has its eigenvalues \j as its diagonal elements and has the 
properties that (i) If S and T are upper triangular matrices with eigenvalues Xsj and Xtj, 
then ST is also a triangular matrix, with diagonal elements {ST)jj = XsjXtj- A corollary 
is that for an arbitrary N x N matrix T 

N 

i=i 

The dimension of Tz,A,Ly,d, or equivalently, TT,A,Ly,d, is [12] 
nz{A,Ly,d)=nT{A,Ly,d) = (^^^^"^^^^ {l^-d) A = sq,tri,hc. (1.20) 

for < d < Lj, and zero otherwise. The property that the numbers ^^(A, Ly, d) are the 
same for all three lattices A = sq, tri, he, as was shown for A = sq, tri in Ref. [12] and for 
A. — he in. [15]. For this reason we shall henceforth usually revert to our previous notation 
in Ref. [12], setting 

nz{A, Ly,d) = nz{Ly,d) , A = sq,tri,hc (1-21) 

The formal quantity nz{0, d) will appear in some determinant formulas below and is given 
by eq. (1.20) as the Kronecker delta, nz{0,d) — Sdfi- (Below we shall consider self-dual 
strips Go of the square lattice, which have different dimensions riziGo, Ly, d); for this case 
we shall include the Go dependence in the notation.) Special cases of nz{Ly, d) that are of 
interest here include 

nz{Ly,Ly) = l (1.22) 

nz{Ly,Ly-l) = 2Ly-l (1.23) 

nz{Ly,Q)^CLy (1.24) 
where here C„ is the Catalan number which occurs in combinatorics and is defined by 
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C„ = -^f^^V (1.25) 
{n + l)\n J 

(No confusion should result from our use of the same symbol Cn to denote the circuit graph 
with n vertices since the meaning will be clear from context) . The first few Catalan numbers 
are Ci = 1, C2 = 2, C3 = 5, C4 = 14, and C5 = 42. 

The full transfer matrix Tx,K,Ly-i ^ — Z,T, has a block structure formally specified by 

Ly 

Tx,A,Ly ^ ®UTx,A,Ly,d X^Z,T (1.26) 

d=0 

where the product Y['^x,A,Ly,d means a set of square blocks, each of dimension, of the 
form Xx,A,Ly,d,j times the identity matrix. The dimension of the total transfer matrix, i.e., 
the total number of eigenvalues Xx,A,Ly,d,j, counting multiplicities, is thus 

Ly Ly 

dim{Tx,A,Ly) = Yl dMTx,A,Ly,d) = E L„ d) X ^ Z,T . (1.27) 

d=0 d=0 

As in our earlier work, we define Nz,A,Ly,x as the total number of distinct eigenvalues of 
Tx,A,Ly, X = Z,T, i.e. the sum of the dimensions of the submatrices Tx,A,Ly,d, modulo the 
multiphcity c^'^l This is [12] 

Nz,A,Ly = NT,A,Ly = E ^^(^2" ^ ( r / ^ = ■ (1-28) 

d=0 ^ ^ 

II. TRANSFER MATRIX METHOD 

There are several equivalent ways to calculate Z{G,q,v) or T{G,x,y) for these lattice 
strip graphs. One is to make iterative use of the deletion-contraction relations obeyed by 
T{G,x,y). This yields a generating function whose denominator directly determines the 
Xx,A,Ly,d,j^s where X = Z,T. This iterative method works also for strip graphs with free 
longitudinal boundary conditions. The transfer matrix method for lattice strips of the square 
and triangular lattices with free longitudinal boundary conditions (and free or cylindrical 
transverse boundary conditions) was explained in detail in Ref. [29]. For the free strip with 
width Ly, the bases of the transfer matrix are all of the possible non-crossing partitions of Ly 
vertices. (For the zero-temperature antiferromagnetic Potts model, the non-nearest-neighbor 
requirement is imposed.) The eigenvalues of the transfer matrix for a free strip are the same 
as the eigenvalues of the transfer matrix Tz,A,Ly,d=o in the degree d — subspace (called "level 
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0" subspace in [28]) for the corresponding cyclic strip, and the dimension of this matrix is 
nz{Ly, 0) = C^y, the Catalan number, as in eq. (1.24). This was shown for the free strips of 
the square and triangular lattices in [29] and for the cyclic/Mobius strips in [12] and extended 
to the honeycomb lattice in [15]. As discussed in Refs. [26]- [28] for chromatic polynomials, 
which we generalize here to the full Potts model partition function, the degree d = 1 subspace 
is given by all of the possible non-crossing partitions with a color assignment to one vertex 
(with possible connections with other vertices), and the multiplicity is g — 1 = c(^). This 
follows because there are q possible ways of making this color assignment, but one of these 
has to be subtracted, since the effect of all the possible color assignments is equivalent to 
the choice of no specific color assignment, which has been taken into account in the level 
subspace. In this derivation and subsequent ones we assume that q is an integer > 4 
to begin with, so that the multiplicities are positive-definite; we then analytically continue 
them downward to apply in the region < q < 4: where c^'^^ can be zero or negative. For the 
next subspace we consider all of the non-crossing partitions with two-color assignments to 
two separated vertices (with possible connections with other vertices) . Now the multiplicity 
can be understood by the sieve formula of [26]- [28]. Since the two assigned colors should be 
different, there are q{q — 1) ways of making these assignments. This includes the q possible 
color assignments for each of the two vertices that have been considered in level 1 and hence 
these must be subtracted. In doing this, the no-color assignment was subtracted twice, and 
one of these has to be added back. Therefore, the multiplicity is 

?(g-l)-2g + l = g2_3g+i = c(2) _ (^2.1) 

By the same method, for the non-crossing partitions with the three-color assignment, the 
multiplicity is 

q{q - if - 2q{q - 1) - q^ + 3q - 1 = q^ - 5q^ + 6q - 1 = c^^' . (2.2) 
The multiplicity for the four-color assignment is 

q{q - 1)3 - 2q{q - if - 2q\q - 1) + 3q{q - 1) + 3q^ - 4q + 1 

= {q - l){q^ - Gq"^ + 9q - 1) = c^^K (2.3) 

and so forth for higher values of d. We show the calculation for these multiplicities pictorially 
for 2 < d < 4 in Fig. 1. In general, the multiphcity of the d-color assignment can be 
computed to be c^^^ given in eq. (1.15). We list graphically all the possible partitions for 
Ly — 2 and Ly — 3 strips in Figs. 2 and 3, respectively, where white circles are the original 
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Ly vertices and each black circle corresponds to a specific color assignment. We remark 
that the connections between the black circles and white circles also obey the non-crossing 
restriction, so it is possible that two non-adjacent black circles represent the same color. 
In the following discussion, we will simply use the names white and black circles with the 
meaning understood. We denote the partitions VLy,d ior 2 < Ly < i as follows: 



7^2,0 = {I; 12} , V2,i = {2; 1; 12} , 7^2,2 = {1, 2} 



(2.4) 



Vs,o = {/; 12; 13; 23; 123} , Vs,! = {3; 2; 1; 12, 3; 12; 13; 23; 23, 1; 123} , 



7^3,2 = {2, 3; I, 3; 1, 2; 12, 3; I, 23} , 7^3,3 = {1, 2, 3} 



(2.5) 



7^4,0 = {/; 12; 13; 14; 23; 24; 34; 12, 34; 14, 23; 123; 124; 134; 234; 1234} , 

^ = {4; 3; 2; i; 12, 4; 12, 3; l2; 13, 4; T3; Ti; 23, 4; 23; 23, 1; 24; 24, 1; 34; 34, 2; 34, 1; 12, 34; 
34, T2; 23, 14; 123, 4; T23; T24; T34; 234; 234, 1; T234} , 
7^4,2 = {3, 4; 2, 4; 1, 4; 2, 3; 1, 3; 1, 2; 12, 3, 4; T2, 4; T2, 3; 13, 4; 23, 4; 23, 1, 4; 1, 23; 1, 24; 2, 34; 



1, 34; 34, 1, 2; 12, 34; 123, 4; 1, 234} , 



7^4,3 = {2, 3, 4; 1, 3, 4; 1, 2, 4; 1, 2, 3; 12, 3, 4; 1, 23, 4; 1, 2, 34} , 



P4,4 = {1,2,3,4} 



(2.6) 
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FIG. 1. Multiplicities of the transfer matrices for 2 < d < 4. 
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FIG. 2. Partitions for the Ly = 2 strip. 
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FIG. 3. Partitions for the Ly = 3 strip. 



Using these methods, we can obtain nz{Ly, d) given in Theorem 4 of [12] as follows. 
Since the maximum number of colors to assign is Ly for a strip with width Ly, it follows 
that nz{Ly,d) = for d > Ly and nz{Ly, Ly) = 1. It is elementary that nz(l,0) = 1. The 
d = partitions of a width- (Lj^ + 1) strip can be obtained by either adding a unconnected 
white circle to the bottom of the d = partitions of a width- L,^ strip or converting the black 
circle of the d = 1 partitions of a width-Lj^ strip into a white circle. That is, nz{Ly + 1,0) = 
nz{Ly,0) + nz{Ly, 1). Finally, the partitions of a width-(Lj^ + 1) strip for 1 < c? < L^^ + 1 
can be obtained in one of the following four ways: (a) adding a pair of connected circles, 
one black and one white, (but not connected to any other vertex) above the highest black 
circle of the d — 1 partitions of a width-Lj^ strip; (b) adding a unconnected white circle 
above the highest black circle of the d partitions of a width-L^^ strip; (c) adding a white 
circle above the highest black circle of the d partitions of a width- L^^ strip and connecting 
these two circles; (d) converting the highest black circle of the d+1 partitions of a width-Lj^ 
strip into a white circle. Now the lowest white circle of the d partitions of a width- (Lj^ -|- 1) 
strip can either connect to a black circle with or without ((c) or (a)) other connections to 
other white circles, or it does not connect to a black circle with or without ((d) or (b)) other 
connections to other white circles. Therefore, (a) to (d) exhaust all the possibilities, and we 
have nz{Ly + l,d) — nz{Ly, d — 1) -\- 2nz{Ly, d) -\- nz{Ly, d -\- 1) lox 1 < d < Ly -\- 1. This 
completes the proof. Indeed, these are precisely the equations that we obtained in a different 
manner in Ref. [12]. Solving these equations as we did in [12], it follows that the dimension 
of the transfer matrix Tz^A,Ly,d is the expression for nz{Ly, d) given in eq. (1.20). 

Similarly, for the zero-temperature antiferromagnetic Potts model, np{A, Ly, d) for 
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A = sq,tri given in Theorem 2 of [12] can be shown as follows. Now the adjacent white 
circles cannot connect to each other. Again, it is obvious that np{Ly,d) = for o? > Ly, 
np{Ly, Ly) = 1, and np(l,0) = 1. The d = partitions of a width-(Lj^ + 1) strip can all 
be obtained by converting the black circle of the d = 1 partitions of a width-Lj^ strip into a 
white circle with the condition that the connection between the black circle and the lowest 
white circle, if there is one, must be deleted. That is, np{Ly + 1,0) = np{Ly, 1). Finally, the 
partitions of a width-(Ly + 1) strip for 1 < d < L + 1 can be obtained in one of the following 
four ways: (a) adding a pair of connected circles, one black and one white, (but not connected 
to any other vertex) above the highest black circle of the d—1 partitions of a width-Lj^ strip; 
(b) adding a unconnected white circle above the highest black circle of the d partitions of 
a width- strip; (c) if the highest black circle of the d + 1 partitions of a width- L^^ strip 
connects to the lowest white circle, deleting this connection, converting this black circle into 
a white circle, and connecting it to the black circle below it; (d) if the highest black circle of 
the d+ 1 partitions of a width- Lj^ strip does not connect to the lowest white circle, converting 
it into a white circle. Now the lowest white circle of the d partitions of a width-(Ly + 1) strip 
can either connect to a black circle with or without ((c) or (a)) other non-adjacent connec- 
tions to other white circles, or it does not connect to a black circle with or without ((d) and 
(b)) other non-adjacent connections to other white circles. Therefore, (a) to (d) exhaust all 
the possibilities, and we have np{Ly + l,d) — np{Ly, d — 1) + np{Ly, d) -\- np{Ly, d+ 1) for 
1 < d < L + 1. This completes the proof. 

The construction of the transfer matrix for each level (= degree) d can be carried out 
by the methods similar to those for d = 0. For cyclic strips, the transfer matrix Tz,A,Ly,d is 
the product of the transverse and longitudinal parts, Hz,K,Ly,d and Vz,K,Ly,d-i which can be 
expressed as 

Ly-l 

Hz,sq,Ly,d = Hz,tri,Ly,d = H ("^ + Ly,d,i,i+l) 

i=l 

[Ly/2] [iLy-l)/2] 

Hz,hc,Ly,d,l — n (-'^ + vJLy,d,2i-l,2i) , Hz,hc,Ly,d,2 — Y\. V J Ly,d,2i,2i+\) 

1=1 1=1 



Vz,sq,Ly,d — Vz,hc,Ly,d — Y\.^^^ ~'~ DLy^(i,i) 

i=l 

Ly-l 

Vz,tr^,Ly,d = n [(^^ + ^Ly,d,)il + vJLy4,,+l)]ivI + DLy,d,Ly) , (2.7) 
i=l 

where [u] denotes the integral part of u, and / is the identity matrix of dimension nz{Ly, d). 
JLy,d,i,i+i is the join operator between vertices i and and Dl^^^i is the detach operator 
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on vertex i [29] . We have 




Tz,tri,Ly,d — yZ,tri,Ly,dHz,tri,Ly,d 



Tz,hc,Ly,d — {Vz,hc,Ly,dHz,hc,Ly,d,2){Vz,hc,Ly,dHz,hc,Ly,d,l) = Tz,hc,Ly,d,2Tz,hc,Ly,d,l ■ 



(2.8) 



Here we shall give results for this method for both cyclic and Mobius strips. 

Consider two sets of Ly vertices and denote them as ii,i2, ■■■''■Ly and ji,j2, ■■■jLy- For cyclic 
strips of the square lattice, the horizontal edges connecting these vertices are (ii, ji), (^2, J2), 

{iLy,jLy)- For Mobius strips, one set of horizontal edges becomes («i,jLy)i (^2ijLy-i)) 
(^Ly) ii)- This corresponds to exchanging the pair of bases that switch to each other when 
the vertices reverse in order, i.e., the set of bases that do not have self- reflection symmetry 
with respect to the center of the Ly vertices. For example, among the partitions for the 
Ly — 2 strip in Fig. 2, the first partition 2 and the second partition I in 7^2,1 niust be 
exchanged under this reflection. The pairs of partitions for the Ly — 3 strip in Fig. 3 are 
(1) the second partition 12 and the fourth partition 23 in Vsfi, (2) the first partition 3 and 
the third partition 1 in 7^31, (3) the fifth partition 12 and the seventh partition 23 in Vs^i, 
(4) the fourth partition 12, 3 and the eighth partition 23, 1 in ^3^1, (5) the first partition 2, 3 
and the third partition 1, 2 in 7^3 2, and (6) the fourth partition 12, 3 and the fifth partition 
1, 23 in 7'3^2- For this specific set of edges of Mobius strips, the pairs of columns of Vz,A,Ly,d 
that correspond to these pairs of partitions should be exchanged, and these matrices will be 
denoted as Vz,A^Ly4- Equivalently, the same pairs of columns of Tz^A,Ly,d should be exchanged, 
and these matrices will be denoted as Tz,K,Ly,d = Vz,A,Ly,dHz,A,Ly,d for A = sq, tri. There 
are two kinds of Mobius strips for the honeycomb lattice. When Ly is even, the number of 
vertices in the horizontal direction is even as for the cyclic strips, i.e., L^ = 2m. When Ly is 
odd, the number of vertices in the horizontal direction is odd, L^ = 2m — 1. Therefore, for 
the honeycomb lattice, we use the definition 



As was discussed in [15] for the crossing-subgraph strips, the square of each eigenvalue of 
TzMcLyxi for odd Ly is an eigenvalue of Tz,hc,Ly,d- 

We have found in [12] the changes of coefficients for the square lattice when the longi- 
tudinal boundary condition is changed from cyclic to Mobius. In the context of our present 
transfer matrix formalism, we can express these changes of coefficients for the square, trian- 
gular and honeycomb lattices as follows: 



Z,hc,Ly,d * Z,hc,Ly,d-'-^ Z,hc^Ly,d,l 



for odd Ly 



rp TT -tr TJ 

-'■ Z,hc,Ly,d * Z,hc,Ly,d-'-'- Z,hc,Ly,d,2* Z,hc,Ly,d-'-'- Z,hc,Ly,d,l 



for even L, 



y 



(2.9) 
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(2.10) 



1< A; < 



2 J 



(2.11) 



^(2fe+l) ^ 



< A; < 



The Potts model partition function for Mobius strips is given by 
Z{A,Lx L,,Mb,q,v) = c^'^Tr[iTz,A,Lyflr~'Tz,A,L„o] 

[iLy-l)/2] 

+ c('^+^)Tr[(Tz,A,L„2.+i)™-^Tz,A,L„2.+i] 

[Ly/2] 

- J2 c('-'^7^r[(rz,A,L„2.)'"-^Tz,A,Z,,,2.] 
d=l 

The corresponding formula for the Tutte polynomial is 

T(A,L X L,,Mb,x,y) = -1- (^c(°Tr[(rT,A,L„o)'"-'7V,A,L„o] 

l{Ly-l)/2] 

+ c('^^'^r^[(rT,A,L„2.+l)"^-^rr,A,L„2d+l] 
d=0 

[Ly/2] 

- J2 c('-'^7^r[(rr,A,L„2.)"^-^T^,A,L„2.]) . 

d=l 



(2.12) 



(2.13) 



(2.14) 



For the square lattice or the honeycomb lattice with Ly even, the eigenvalues of Tz,K,Ly,d 
are the same as those of Tz^A,Ly,d except for possible changes of signs. The number of 
eigenvalues with sign changes is equal to the number of column-exchanges from T^^A.Ly.d to 
Tz,A,Ly,d- Denote the number of eigenvalues that are the same for Tz,sq,Ly,d and Tz,sq„Ly,d as 
nz{sq, Ly, d, +), and the number of eigenvalues with different signs as nz{sq, Ly, d, —). It is 
clear that 



nz{sq, Ly, d, +) + nz{sq, Ly, d, -) = nz{sq, Ly, d) . 



(2.15) 



Define 



Anz{sq, Ly, d) = nz{sq, Ly, d, +) - nz{sq, Ly, d, -) 



(2.16) 
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which gives the number of partitions which have self-reflection symmetry. For example, 
among the partitions for the Ly = 2 strip in Fig. 2, the partitions / and 12 in 7^2,0) the third 
partition 12 in 7^2,1) and the partition 1,2 in 7^2,2 have self-reflection symmetry. Among the 
partitions for the Ly = 3 strip in Fig. 3, it includes the first partition /, the third partition 
13 and the fifth partition 123 in V^^, the second partition 2, the sixth partition 13 and the 



ninth partition 123 in Vs^i, the second partition 1,3 in 7^3,2, and the partition in 7^3,3. We 
list Anz{sq, Ly,d) for 1 < L^, < 10 in Table I in the appendix. The total number of these 
partitions for each Ly will be denoted as ANz^Ly = '2^'"- The relations between Anz{sq, Ly, d) 
are 

Anz{sq, 2n, 0) = 2Anz{sq, 2n - 1, 0) , for < n 
Anz{sq, 2n, 2m — 1) = Anz{sq, 2n, 2m) 

= Anz{sq, 2n — 1, 2m — 1) + Anz{sq, 2n — 1, 2m) , for 1 < m < n 
Anz{sq, 2n+ 1, 2m) = Anz{sq, 2n+ 1, 2m + 1) 

^ Anz{sq,2n,2m) + Anz{sq,2n,2m + 1) , for < m < n . (2.17) 

We also list nz{sq, Ly, d, +) and nz{sq, Ly, d, — ) for 2 < Ly < 10 in Table II in the appendix. 
Notice that nz{sq,Ly,{),+) is the number of Xz,sq,FF,Ly proved in [30]. Recall the numbers 
of \z,sg,Ly,j for the Mobius strips of the square lattice with coefficients ic*-^-*, defined as 
nz,Mb{Ly, d, ±) = nz,Mb{sq, Ly, d, ±), have been given in [12]. With the eqs. (2.10) to (2.12), 
the relations between nz{sq, Ly, d, ±) and nz,Mb{sq, Ly, d, ±) are 

nz,Mb{sq, Ly, 0, ±) = nz{sq, Ly, 0, ±) + nz{sq, Ly, 2, ^) 

nz,Mb{sq, Ly, k, ±) = nz{sq, Ly, 2k -1,±) + nz{sq, Ly, 2k + 2,^) , 



for l<k< 



Ly+1 



2 J- ^^-''^ 
For a strip graph Gs of a lattice A with given boundary conditions, following our earlier 
notation [11] we denote the sum of coefficients (generalized multiplicities) Cq^ of the Xz,Gs,j 
as Cz,Gs- This sum is independent of the length m of the strip. For a cyclic strip graph this 
is given by [12] 

Ly 

Cz,A,cyc.,Ly = c^'^^nziLy, d) = q^y for A = sq, tri, he . (2.19) 

d=0 
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For the Mobius strip of the square lattice or the honeycomb lattice with Ly even, the sign 
changes of the eigenvalues of Tz^A,Ly,d can be considered as the sign changes of the coefficients. 
For these cases, the sum of coefficients is given by Theorem 8 of Ref. [12]: 



n — ST^ — A IT A\ (d) _ ) 0.^' even L„ 

<^Z,(sq,hc),Ly,Mb = 2^ (^Z,Ly,Mb,j - /\nz,Mb[-l^y, a)C - < (Ly+l)/2 f^j. ^ 

(2.20) 



where 

dma.-[t^,, for even 

Y^y^ for odd ^ ' 

In previous work we have given results for the determinants for various strip graphs Gg 
(e.g., [13]) 

detrz(G.)= n i>^z,G.,or'°-^ (2.22) 

and 

detrp(a)= n (^P,G.J'"'°^'^^ ■ (2.23) 

and we shall extend these results to arbitrary width here. In the present context, these can 
be written, for cyclic strips, as 

Ly 

det{Tz,K,Ly) = \{[det{Tz,K,Ly4)f' (2.24) 

det{TT,A,Ly) = Ylldet{TT,A,Ly,a)f' . (2.25) 

d=0 

From eq. (1.12) it follows that [11,15] 

>^Z,A,Ly,d,3 = V^^" >^T,A,Ly,d,j (2.26) 

and 

Tz,A,Ly,d = vP'^yTT,A,Ly,d (2.27) 

where 

_ f 1 if A = or tri or Gd /q qqn 

^"12 iiA = hc ^^^^ 
Note that the factor of {x — 1) in eq. (1.12) cancels the factor l/{x — 1) in eq. (1.17). 
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III. PROPERTIES OF TRANSFER MATRICES AT SPECIAL VALUES OF 

PARAMETERS 

In this section we derive some properties of the transfer matrices Tz,K,Ly,d at special values 
of q and and, correspondingly, TT,K,Ly,d at special values of x and y. 

A. ^ = 

Prom (1.1) or (1.8) it follows that for any graph G, the Potts model partition function 
Z{G, q, v) satisfies 

Z(G,g,0) = g-(^) . (3.1) 

Since this holds for arbitrary values of g, in the context of the lattice strips considered here, 
it implies 

{Tz,K,Ly,d)v=o^^ for l<d<Ly (3.2) 

i.e. these are zero matrices. Secondly, restricting to cyclic strips for simplicity, and using the 
basic results n = LyL^ = Lym for A — sq,tri and n — 2Lym for A = /ic, eq. (3.1) implies 
that 

rp^rrp ^m^ _ / q^""' for A = sq, tri „x 

J r[{lz,K,cyc.,Ly,q,v) \v=0 - | ^2Ly7n K = hc ^ 

With our explicit calculations, we find that 

(Tz,A,L„o)ifc = for ^ = 0> and j>2 (3.4) 

i.e., all rows of these matrices except the first vanish, and 

(Tz,A,L„o)ii = g^"-^ for ^; = (3.5) 

where p was given in eq. (2.28). Por this v — Q case, since all of the rows except the first 
are zero, the elements {Tz,K,Ly,Q)ik for A; > 2 do not enter into Tr[{Tz,K,Lyfi)^\-, which just 
reduces to the m'th power of the (1,1) element: 

TT[{Tz,^,:^^,or] = [{Tz,K,Ly,o)iir ■ (3-6) 

Corresponding to this, all of the eigenvalues Xz,K,Ly,d,j vanish except for one, which is equal to 
{Tz,!^,Ly,o) 11 in eq. (3.5). As will be seen, this is reflected in the property that det{Tz,K,Ly,d) 
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has a nonzero power of v as a factor for Ly > 2 for all of the lattice-A strips considered 
here. The restriction > 2 is made because the strips of the triangular and honeycomb 
lattice are well-defined without degenerating for Ly > 2 and, in the case of the square lattice, 
for the case Ly = 1, both of the transfer matrices reduce to scalars, Tz,sq,ifi = q + v and 
Tz,sq,i,i = V, the former of which is, in general, nonzero at v = 0. Note that the condition 
V = is equivalent to the Tutte variable condition y = I. 

B. q = 

Another fundamental relation that follows, e.g., by setting q ~ in eq. (1.8), is 

Z(G,0,v)^0. (3.7) 
Now the coefficients c^'^^ evaluated at q = satisfy [12] 

c('^) = (-l)<^ for g = . (3.8) 
Hence, in terms of transfer matrices, we derive the sum rule 

E Tr[{Tz,A,L,,,r] - Yl Tr[{Tz,A,L,,,r] =0 for g = . (3.9) 

0<d<Ly, d even ^<d<Ly, d odd 

This is similar to a sum rule that wc obtained in Rcf. [21]. Since it applies for arbitrary m, 
it implies that there must be a pairwise cancellation between various eigenvalues in different 
degree-d subspaces, which, in turn, implies that at g = there are equalities between these 
eigenvalues. For example, for Ly = 1 and q = 0, A^^^g^i^o = q+'v becomes equal to Xz,sg,i,i — v. 
For the Ly — 2 strips of all of the A = sq, tri, he lattices, g = 0, two of the eigenvalues of 
2z,A,2,i become equal to the eigenvalues of Tz,k,2,o-, while the third eigenvalue of Tz^k,2,\ 
becomes equal to \z,k,2,2 — f and so forth for higher Ly. Note that setting g = does not, 
in general, lead to any vanishing eigenvalues for the Tz,K,Ly,d and hence our formulas below 
for Tz,A,Ly,d do not contain overall factors of g. 

In terms of Tutte polynomial parameters, the value g = implies x — 1 (unless v — Q). In 
contrast to the vanishing of Z{G, g, i^) at g = 0, the Tutte polynomial T{G, 1, y) is nonzero for 
general y. This different behavior can be traced to the feature that in eq. (1.12), Z{G, g, v) 
is proportional to T{G,x,y) multiplied by the factor [x — 1)'^('^) — (^x — 1), so at x = 1, 
Z{G, 0, -u) = even if T{G, 1, y) 0. 
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C. q = l 



Evaluating eq. (1.1) for q = I, one sees that the Kronecker deha functions Sa^aj = 1 for 
aU pairs of adjacent vertices consequently, 



Z{G, 1, v) = e^(^)^ = (1 + v)<^^ . (3.10) 



Now [12] 



r 1 if = mod 3 
If q^l then c^*^) = <^ if d = 1 mod 3 (3.11) 

I -1 if d = 2 mod 3 

Hence, in terms of transfer matrices, we derive the sum rule for the present cyclic lattice 
strips G — A, Ly X Lx, eye. 

E Tr[{Tz,A,L,,ar] - E Tr[{Tz,A,L,,ar] = (1 + vT^''^ for g = 1 , 

0<d<Ly,d=0 mod 3 2<d<Ly,d=2 mod 3 

(3.12) 

Here the number of edges e{G) for each type of cychc strip is 

{{2Ly - l)m if A = 
'^Ll'-t lilt' ("3) 
2Lym if A = 

where m is given in terms of by eq. (1.4). As before, the sum rule (3.12) implies relations 
between the eigenvalues of the various transfer matrices Tz,K,Ly,d- 

Y). V = —1, i.e., y = 

As discussed above, the special value v = —1, i.e., y = 0, corresponds to the zero- 
temperature Potts antiferromagnet, and in this case the Potts model partition function 
reduces to the chromatic polynomial as indicated in eq. (1.7). In Refs. [12,17,15] we have 
determined how the dimensions nz{Ly,d) for the square, triangular, and honeycomb lattice 
strips and nziGo, Ly,d) for the self-dual strip of the square-lattice reduce to dimensions 
np{Ly,d) and np{G]:), Ly,d), respectively. In all cases except for the lowest width Ly = 1 
for the square lattice strip, in each degree-d subspace for < ci < — 1 for A = sq, tri, he 
and 1 < d < for A = Go, some eigenvalues vanish so that np{A, Ly,d) < nz{A, Ly,d) 
and hence det{Tz,K,Ly,d) — 0. This property is reflected in the powers of {v + 1) and y that 
appear, respectively, in our formulas below for det{Tz,A,Ly,d) and det{TT,K,Ly,d)- 
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E. g = —V, i.e., x = 



For the graph G = G{V,E), setting x = and y = 1 — g in the Tutte polynomial 
T{G, X, y) yields the flow polynomial F{G, q), which counts the number of nowhere-0 g-fiows 
(without sinks or sources) that there are on G [10]: 

F{G, q) = {-1)<G)-^(G)+^T{G, 0,1 - q) . (3.14) 

In Ref. [38] we showed that the flow polynomial for cyclic lattice strips has the form 

Ly np(A,Ly,d) 

F{A, Ly X L., eye, g) = c^"^ ^^^AMjT (3-15) 

d=0 3=1 

or equivalently, 

Ly nF{A.,Ly,d) 

F{A, Ly X L,, eye, q)^Yl ^^'^ E ^^[(^^,A,L„dj)"^] ■ (3.16) 

d=0 j=l 

We determined the dimensions ^^(A, Ly, d) and found that for all Ly > 1, npi-A, Ly, d) < 
nz{A.,Ly,d). Thus, when one sets x = in the Tutte polynomial, or equivalently, q — —v 
in the Potts model partition function, some of the eigenvalues in each degree-d subspace 
vanish, and hence det{Tz,K,Ly,d) and det{TT,K,Ly,d) vanish. This is reflected in the powers 
of {q + v) and x that appear, respectively, in our formulas below for det{Tz,K,Ly,d) a-nd 
det{TT^A^Ly,d)- In passing, we note that for planar graphs Gpi, because of the duality relation 
(1.13), the chromatic polynomial on G is closely related to the flow polynomial on G*: 
P{G,q) = qF{G\q). 

These results concerning evaluations of the Potts model partition function on these lattice 
strips for special values of q,v, equivalent to evaluations of the Tutte polynomial for special 
values of x and y, help to explain the types of factors that appear in the simple formulas 
that we obtain below for the determinants det{Tz^A,Ly,d) and det{TT^A,Ly,d)- 

IV. GENERAL RESULTS FOR CYCLIC STRIPS OF THE SQUARE LATTICE 

In this section and the subsequent ones we present general results that we have obtained, 
valid for arbitrarily large strip width Ly (as well as arbitrarily great length) for transfer 
matrices and their properties. We begin with the strip of the square lattice. 
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A. Determinants 



We find 

det{TT,s,,Ly4) = (4.1) 

where nz{Ly,d) was given in eq. (1.20). This apphes for all d, i.e., < d < Ly since [12] 
nz{Ly — l,d) = ioi d > Ly — 1. This is equivalent to the somewhat more complicated 
expression for det{Tz,sq,Ly,d)- 



det{Tz,s,,Ly,d) = {v^^r^^^^'^ [(I + + vfy-'\ " . (4.2) 

These determinant formulas can be explained as follows. By the arrangement of the partitions 
as shown in Figs. 2 and 3, the matrix JLy,d,i,i+i has the lower triangular form and the matrix 
Di^^ ,i,i has the upper triangular form. From the definition of the transfer matrix in eq. 
(2.8), the determinant of Tz,sq,Ly,d is the product of the diagonal elements of Hz,sq,Ly,d and 
Vz,sq,Ly,d- The diagonal elements of Hz,sq,Ly,d has the form (1 + vY , where r is the number 
of edges in the corresponding partition, and the diagonal elements of Vz.sq.Ly.d has the form 
v^y(l + q/vY ^ where s is the number of vertices which does not connect to any other vertex 
in the corresponding partition. Therefore, the power of (1 + v) in eq. (4.2) is the sum of 
the number of edges of all the (L^,, rf) -partitions. Let us compare the (L^,, rf) -partitions and 
{Ly — 1, d) partitions. Since each edge of a {Ly^ (i)-partition corresponds to converting a vertex 
of a {Ly—1, (i)-partition into a edge, the power of (1+";-^) is {Ly — l)nz{Ly—l, d). It is clear that 
the power of v is Lynz{Ly, d). The power of (1 + q/v) in eq. (4.2) is the sum of the number 
of unconnected vertices of all the {Ly, (i)-partitions. Compare the {Ly, d) -partitions and 
{Ly — 1, d) partitions again. Since each unconnected vertex of a {Ly, (i)-partition corresponds 
to adding an unconnected vertex to a {Ly — 1, d)-partition with Ly possible ways (between 
every two adjacent vertices, above the highest vertex, and below the lowest vertex), the 
power of (1 -|- q/v) is Lynz{Ly — l,d). 

Next, taking into account that the generalized multiplicity of each Xz,\,Ly,d is c^'^\ we 
have, for the total determinant 

Ly 



nz(Ly-l,d) 



det{Tz,sq,Ly) = J[[det{Tz,sq,Ly,d)X 

d=0 

Ly 

^ Yl^^y ~ l)^'']"^(^^''^)'=^'"[a;^^2/^^"^]''^(^^~^''^)''^''^ 

d=0 

^ [(y _ i^Ly-^Edlo^z(I^y'd)<^^''^[x^yyI^y-^lEdlo''z{Ly-l,d)c(''^ _ (4 3) 
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Using eq. (2.19) together with eq. (4.2) from Ref. [12], viz., nz{Ly,d) ~ for d > Ly, so 
that Y!d=Q'^z{Ly - 1, c?)c('^) = Sdlo^ nz{Ly - 1, d)c^'^\ we have, finally, 

det{Tz,s,,L,) = [{y - 1)^1'^'^ [x'^yy'^y-Y'''' ■ (4-4) 

This agrees with the conjecture given as eq. (2.42) of our earlier Ref. [13] for the determinant 
of the transfer matrix of the cyclic strip of the square lattice with arbitrary width Ly. 



B. Traces 



For the trace of the total transfer matrix, taking account of the fact that each of the 
Xx,A,d,j, X = Z,T, has multiplicity c^'^\ we have 



Ly nz{Ly,d) ^ Ly 

d=0 7=1 d=0 



d=0 



(4.5) 



Note that Yjd=o'^'^^^'^''"i'^T,sq,Ly,d) contains a factor of {x — 1) which cancels the prefactor 
l/{x — 1). Equivalently, we find 



Tr{Tz,sq,Ly) = J2 ^^'^ E ^z,sq,Ly,d,j = E c^'^Tr{Tz,sq,Ly,d) 

d=0 3=1 d=0 



(4.6) 



We derive this as follows. The trace here is the m — 1 case in eq. (1.17) or (1.6), which 
corresponds to a L^^-vertex tree with a loop attached to each vertex, as shown in Fig. 4. The 
corresponding Tutte polynomial is x^y^^y^y. In principle, we can also consider the m = 1 
case in eq. (2.13) or (2.14) for the Mobius strips, but the result is not as simple as that listed 
here. 



2^71 — I9 ^Lii — 1 



Ly-1 



Jy ±Jy 

FIG. 4. m = 1 graph for the cyclic square lattice. 
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C. Eigenvalue for d = Ly for A = sq,tri, he 



It was shown earlier [11,39] that the A's are the same for a given lattice strip with cyclic, 
as compared with Mobius, boundary conditions (and it was shown that the A's for a strip 
with Klein bottle boundary conditions are a subset of the A's for the same strip with torus 
boundary conditions). Prom eq. (1-20) one knows that there is only one A for degree d = Ly, 
which we denote \z,K,Ly,Ly- That is, for this value of d, the transfer matrix reduces to 
1x1, i.e. a scalar. We found that for a cyclic or Mobius strip of the square, triangular, or 
honeycomb lattice with width Ly, 

^T,K,Ly,Ly = 1 ■ (4.7) 

(Indeed, we have found that this holds more generally for other strips of regular lattices 
(e.g., [19]), but we shall not pursue this here.) Note the universality; i.e., this is the same 
for all of these lattices. In terms of Potts model variables the results do not maintain the 
full universality: 

>^z,K,Ly,Ly = for A = sq, tri (4.8) 
Xz,A,Ly,Ly ^ v^'^y for A^hc. (4.9) 

D. Transfer Matrix for d = Ly — 1, A = sq 

Prom eq. (1.20) it follows that for A = sq,tri or he, the number of Xx,Ly,d,j, X — Z or 
X — T , ior d — Ly — 1 is nz{Ly, Ly — 1) = 2Ly — 1, i.e. the transfer matrix in this subspace, 
Tx,A,d: X — Z or X — T, is a, (square) {2Ly — 1) -dimensional matrix. We recall that for 
Ly — 1, TT,sq,i,Q — X. For Ly > 2 we find the following general formula, which we write for 



the Tutte polynomial, since it has a somewhat simpler form. 

{TT,sq,Ly,Ly-i)j,j = ^+x for j = l and j = Ly (4.10) 

{TT,s<i,Ly,Ly-i)j,j = 2 + X for > 3 and 2 < J < Lj, - 1 (4.11) 

{TT,Sq,Ly,Ly-lh,j + l = {TT,Sq,Ly,Ly-lh + l,j = 1 fOf 1 < J < - 1 (4.12) 

iTT,sq,Ly,Ly-i)j,j^y for L + l<j<2Ly-l (4.13) 
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{TT,sq,Ly,Ly-l)j,Ly+j - {TT,sq,Ly,Ly-l)j+l,Ly+j - for l<j<Ly-l (4.14) 

{TT,sq,Ly,Ly-l)Ly+j,j = {TT,sq,Ly,Ly-l)Ly+j,j+l = 2/ - 1 for l<J<Ly-l (4.15) 

with all other elements equal to zero. The elements of Tz^sq,Ly-i are given by the relation 
(2.27). Thus, Tx,sq,Ly,Ly-i for X — T,Z consists of four submatrices: 

1. an upper left square submatrix with indices i,j in the ranges 1 < i, j < Ly and nonzero 
elements given by eqs. (4. 10)- (4. 12) 

2. a lower right square submatrix with indices in the ranges Ly + 1 < i, j < 2Ly — 1 and 
nonzero elements given by eqs. (4.13) 

3. an upper right rectangular submatrix with nonzero elements given by eq. (4.14) 

4. a lower left rectangular submatrix with nonzero elements given by eq. (4.15). 

For general x and y, TT,sq,Ly,Ly-i has rank equal to its dimension, 2Ly — 1. In the special 
case y = which yields the chromatic polynomial, the rank is reduced to Ly for general x 
(and may be reduced further for particular x). In the special case x = which yields the 
flow polynomial, the rank is reduced to — 1 (and may be reduced further for particular 



We illustrate these general formulas for the cases Ly = 2,3,4. For this purpose we 
introduce the abbreviations 



y)- 



y 



(4.16) 



y 



- 1 



and use v — y — 1 ss before. 




(4.17) 



,sg,3,2 — 



I Xi I \ 

I X2 I Td Td 

1 Xi Td 

V V y 



(4.18) 



\ V V y J 
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/ xi 1 rd \ 

1 ,T2 1 ra Td 

I X2 I ra Td 

1 xi 

V y 

w t; 1/ 



(4.19) 



\QQvvOOy/ 



Thus, in general, the upper left-hand submatrix has a main diagonal with end entries equal 
to Xi and interior entries equal to X2- Adjacent to this main diagonal are two diagonals 
whose entries are 1, and the rest of the submatrix is comprised of triangular regions filled 
with O's. The upper right-hand submatrix has a band of two diagonals whose entries are 
together with triangular regions filled with O's. The lower left-hand submatrix has a band 
of two diagonals whose entries are v together with triangular regions filled with O's. And 
finally, in the right-hand lower submatrix the entries on the main diagonal are equal to y 
and the rest of this submatrix is made up of triangular regions of O's. For the lowest values 
Ly = 1,2, some of these parts, such as the triangular regions of zeros, arc not present. 

Using eq. (2.27), it is straightforward to obtain the Tz,sq,Ly,Ly-i from these TT,sq,Ly,Ly-i 
matrices. For example. 



Having determined the general form of TT,sq,Ly,d for d = Ly — l,we calculate its eigenvalues 
^T,sq,Ly,Ly-i,j- We find that these consist of one that is the same independent of Ly, namely 




(4.20) 



'T,sq,Ly,Ly-l,l 



(4.21) 



together with 2{Ly — 1) quadratic roots. It is convenient to label these as {j, ±) with 
2 < J < Ly. We find 




where 




(4.23) 



Note that the product {XT,sq,Ly,Ly~i,j,+){>^T,sq,Ly,Ly-i,j-) is independent of j: 



{XT,sq,Ly,Ly-l,j,+){XT,sq,Ly,Ly-l,j,-) ^Xy , for 2 < j < Ly . 



(4.24) 
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This generalizes to the Tutte polynomial our determination of the Ly — 1 terms Xp,sq,Ly,Ly-i,j 
for the chromatic polynomial in eqs. (7.1.2) and (7.1.3) of Ref. [21], and asq,Ly,j was given in 
eq. (7.1.3) of that paper. Note that in the special case y = (w = —1) in which the Tutte 
polynomial or Potts model partition function reduces to the chromatic polynomial, one of 
each of the Ly — 1 pairs of quadratic roots vanishes and the other becomes x + bgq^Lyj for 
the Tutte polynomial, or equivalently, {—l)^^{asq,Ly,j — q) for the Potts model, in agreement 
with eq. (7.1.2) of Ref. [21]. The term {-lfHa.sq,Ly,i - q) = (-1)^^(1 - q) in eq. (7.1.2) of 
Ref. [21] corresponds to the v = —1 special case of eq. (4.21). 

As corollaries of our general result for TT^sq,Ly,Ly-i we calculate the trace and determinant. 
For this purpose, we note that 

E b.,.L,., = 4 i: cos^ (- (■^. + l-J> ) = 2(L„ - 1) (4.25) 

1=2 j=2 y 



det{TT,sq,Ly,Ly-i) = x^-l/^--' (4.26) 
which is a special case of eq. (4.1), and 

TT{TT,sq,Ly,Ly-i) = X + {Ly - l){2 + X + y) = {Ly - 1)(2 + y) + Lyx . (4.27) 
In terms of Potts model quantities these results are 

\z,sq,Ly,Ly-i,i^v'^y-\v + q) (4.28) 



A 



Z,sq,Ly,Ly-l,j± 



q + v{v + bsq,Ly,j + 2)±J[q + v{v + 6,,,^,,, + 2)]^ - Av{v + q){v + 1) 



for 2<j<Ly 



(4.29) 



det{Tz,s 



2Ly(Ly — l) , 



v + qf^{v + iy 



(4.30) 



and 



Tr{Tz ,sq,Ly,Ly 



yLy- 



{4Ly - 3)v + {Ly - l)v^ + Lyq 



(4.31) 
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V. GENERAL RESULTS FOR CYCLIC STRIPS OF THE TRIANGULAR 

LATTICE 



A. Determinants 



We find 



det{TT,tri,L,,a) = (a;^^y2(L,-l))n.(L,-l,.) 



Equivalently, 



-1) 



nziLy — l,d) 



(5.1) 



(5.2) 



By the same argument as before, the powers of v and (l + q/v) are equal to the corresponding 
powers for the square lattice case. Comparing Vz^sq,Ly,d and Vz^tri,Ly,d in eq. (2.7), one sees 
that a set of (J + vJLy,d,i,i+i) has been included for the triangular lattice, so that the power 
of (1 + f ) becomes twice of the corresponding power for the square lattice. 

Taking into account that the multiplicity of each \x,A,Ly,d,j, X = Z,T is c^'^\ it follows 
that the total determinant is 



^Lyy2{Ly-l)y^y- 



det{TT,tri,Ly) = Yl[det{TT,tri,Ly,dW'^ 
d=0 

Equivalently, 

det{Tz,tr^,Ly) ^ Ylidet{Tz,tr^,Ly,d)f' = {v'^^f' [(l + ^)^'^(1 + V 



(5.3) 



2{Ly-l] 



(5.4) 



d=0 



B. Traces 

For the total trace, taking account of the fact that each of the Xx,A,d,j has multiplicity 
we have 



Tr{TT, 



tri.Ly 



Ly nz(Ly,d) 



d=0 j=l 



tri,Ly,d,j 



X 



^X;c('^)Tr(T^,„,^,,,) 

d=0 



(5.5) 



Equivalently, 
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Tr{Tz,tri,Ly) = J2 ^^'^^ Yl ^Z,sq,Ly,d,j = Yl C^'^^Tr{Tz,tri,Ly,d) 
d=0 3=1 d=0 

= q{q + 2v + v^)^^~^{v + 1)^^ . 



(5.6) 



Again the trace here is related to the m — 1 case in eq. (1.17) or (1.6), which corresponds to 
a Lj^-vertex tree with each edge doubled and with a loop attached to each vertex as shown 
in Fig. 5. Therefore, the Tutte polynomial of this graph is {x + y)^y~^y^^. 



Ly-1 



Ly-1 

Li, 



Ly-U 

L,, 



FIG. 5. m = 1 graph for the cyclic triangular lattice. 



C. Transfer Matrix for d = Ly — 1, A = tri 

As before, from eq. (1.20), we know that the dimension of this transfer matrix is 2Ly — 1. 
The first nontrivial case is = 2. For Ly > 2 we find the following general formula, which 



we write for the Tutte polynomial, since it has a somewhat simpler form. 

{TT,tri,Ly,Ly-l) 1,1 = 2 + x + y (5.7) 

{TT,tri,Ly,Ly-ih,j^S + x + y for 2<j<Ly-l (5.8) 

{TT,tri,Ly,Ly-l)Ly,Ly = 1 + ^ (5.9) 

{TT,tri,Ly,Ly-i)i+j,j ^ 1 for l<j<Ly-l (5.10) 
{TT,tri,Ly,Ly-i)j,k = ^ + x + 2y for 1 < j < Ly - 2 and j + l<k<Ly-l (5.11) 

{TT,tri,Ly,Ly-i)j,Ly = l + x + y for l<j<Ly-l (5.12) 

{TT,tri,Ly,Ly-l)l+J,Ly+J ^ for l<j<Ly-l (5.13) 
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{TT,tr^,Ly,Ly-l)J,Ly+k = ^ for 1 < j < Ly - 1 and 3<k<Ly-l (5.14) 

{TT,tri,Ly,Ly-i)j,k^{y-l){2 + y) for l<j<Ly-2 and j + l<k<Ly-l (5.15) 

{TT,tri,Ly,Ly-i)Ly+j,i+k^ {y-'^){^ + x + 2y) for l<j<Lj/-2 and j < k < Ly - 1 

(5.16) 

(TT,i..,L„L,-i)L„, = (|/-l)(l + a; + y) for 1 < j < - 1 (5.17) 

{TT,tri,Ly,Ly-i)j,Ly+k = y{'^ + y) for 1 < j < A; < - 1 (5.18) 

with all other elements equal to zero. The elements of Tz,tri,Ly,Ly-i are given by the relation 
(2.27). Thus, Tx,tri,Ly,Ly-i can again be usefully viewed as consisting of various submatrices. 

For general x and y, TT,tri,Ly,Ly-i has rank equal to its dimension, 2Ly — 1. In the special 
case y = (chromatic polynomial), the rank is reduced to Ly for general x (and may be 
reduced further for particular x). In the special case x — Q (flow polynomial), the rank is 
reduced to 2{Ly — 1) (and may be reduced further for particular y). 

We illustrate these general formulas with some explicit examples for Ly — 2,3,4. For 
compactness of notation, we use the abbreviations 

tii=2 + x + y, lud = ^ + x + y, lu = 3 + x + 2y, (5.19) 
h^iy- l)lu ={y-l){3 + x + 2y), U^{y- l){y + 2), (5.20) 
c„ = l + x + y, Q = (y-l)c„ = (y-l)(l + x + y), (5.21) 

r„ = ^^^±il, Ti = y{y+l) (5.22) 

y- 1 

(with xi — l-\-x and ra = y/{y — I) as above). Then 

til Cu 

TT,tri,2,i =1 I ra \ (5.23) 

Id ci ri 
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T,tri,3,2 



( tx\ la Cu Tu Tu \ 

1 Ld C^a Tu 

1 x-1 Td 

Id k ci ri n 

V Id ci Q n J 



(5.24) 



T,tri,A,Z 



^ ^11 "^u '^u ^ 

1 I'ud I'u C'u ^d 

1 lud Cu Td r„ 

1 ,Ti rd 

Id k k ci ri ri ri 

Id h ci ri ri 

V Id ci ri J 



(5.25) 



It is straightforward to obtain the Tz^tri,Ly,Ly-i from these TT^tri,Ly,Ly-i matrices; for example, 



v{q + 'iv + v'^) v{q + 3v + v^) v{2 + v){l + v) 



Z,tri,2,l 



v{l + v) 



(5.26) 



v^{3 + v) v^{q + 3v + v'^) v'^{2 + v){l + v) 



VI. GENERAL RESULTS FOR CYCLIC STRIPS OF THE HONEYCOMB 

LATTICE 

A. Determinants 

We find 

det{TT,Hc,Ly,d) = {x^^yy^y-'f^^^y-^''') . (6.1) 

Equivalently, 

,r/ n\2Ly -inziLy — l,d) 

det(Tz,Hc,Ly,d) = (v'^yr-^^y'^^[(l + ^) (1 + vfy-'\ . (6.2) 



This can be understood as follows: by an argument similar to that given before, the power 
of (1 + v) is the same as for the square lattice case. Comparing Tz^sq,Ly,d and Tz^hc,Ly,d in 
eq. (2.8), one sees that Vz.hcLyj = yz,sq,Ly,d has been multiplied twice for the honeycomb 
lattice, so that the powers of v and (1 + q/v) become twice of the corresponding powers for 
the square lattice. 

Taking into account that the multiplicity of each Xx,A,Ly,d,j: X = Z,T is c^'^\ it follows 
that the total determinant for the he lattice is 
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det{TTM,L,) ^ Y[[det{TT,Hc,L,,ci)f ' = (x^^^y^-^)^'^" . (6.3) 



d=0 

Equivalently, 

det{Tz,Hc,Lj ^ Y[[det{Tz,Hc,L,,,)f ' = {v^'^^f'' [(l + f + ^)^-f ' ■ (6.4) 

Summarizing the connections between the determinants of the transfer matrices for the 
three lattice strips, det{Tx,tri,Ly,d) is related to det{TT^sq,Ly,d) by the replacement y ^ y"^ 
(holding X fixed), while det{TT^hc,Ly,d) is related to det{TT,sq,Ly,d) by the replacement x ^ 
(holding y fixed). This, together with the fact that nz{Ly,d) is the same for all of these 
three lattices means that the total determinants det{TT^tri,Ly) and det{TT^hc,Ly) are related to 
det{TT^sq,Ly) by the same respective replacements. Correspondingly, det{Tz^tri,Ly,d) is related 
to det{Tz^sq,Ly,d) by the replacement of (1 + v) by (1 + vY and det(Tz^hc.Ly.d) is related to 
det{Tz,sq,Ly,d) by the replacements of the respective factors v by v"^ (cf. eq. (2.28)) and (1 + ^) 

by (1 + ^)< 

B. Traces 

For the total trace, taking account of the fact that each of the \x,K,d,j has multiphcity 
cW, we have 

^ Ly nz{Ly,d) ^ Ly 

Tr{TT,nc,Ly) = ^ E ^^'^ E >^T,hc,Ly,dJ = ^ E <^^'^Tr{TT,nc,Ly,d) 

d=0 j=l d=0 

^x^^-^{x + yf^ . (6.5) 

Equivalently, 

Ly nz{Ly,d) Ly 

Tr{Tz,hc,Ly) =J2^^'^^ E >'z,hc,Ly,d,j^^c^'^^Tr{Tz,hc,Ly,d) 

d=0 j=l d=0 

= q{q + v)^y-^{q + 2v + v^)^^ (6.6) 

The trace here is related to the m — 1 case in eq. (1.17) or (1.6), which corresponds to 
a 2Lj,-vertex tree with each odd edge doubled as shown in Fig. 6. Therefore, the Tutte 
polynomial of this graph is x^y~^{x -\-y)^y. 
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FIG. 6. m = 1 graph for the cyclic honeycomb lattice with even Ly 



C. Transfer Matrix for d = Ly — 1, A = he 

As before, from eq. (1.20), we know that the dimension of this transfer matrix is 2Ly — 1. 
The first nontrivial case is Ly = 2. Recall in eq. (2.8) the transfer matrix for the honeycomb 
lattice, Tz,hc,Ly,Ly-i, is the product of Tz,hc,Ly,Ly-i,i and TzMc,Ly,Ly-i,2- For Ly > 2 we find 
the following general formula, which we write for the Tuttc polynomial. 



{TT,hc,Ly,Ly-i,i)j,j = 1 + X for l<j<Ly-l 



(6.7) 



{TT,hC,Ly,Ly-l,l)Ly,Ly , ^ 



1 + X for Ly even 



for Ly odd 



(6.8) 



{TT,hc,Ly,Ly-l,l)2j-l,2j = {TT,hc,Ly,Ly-l,l)2j,2j-l = 1 for 1 < J < [Ly/2] (6.9) 



iTT,hC,Ly,Ly-l,\)Ly+2j-l,Ly+2j-l - V ^OV 1 <j < [Ly/2] 



(6.10) 



{TT,hC,Ly,Ly-l,l)Ly+2j,Ly+2j = 1 for 1 < J < [{Ly " l)/2] 



(6.11) 



y 

{TT,hc,Ly,Ly-l,l)2j-l,Ly+2j-l = {TT,hc,Ly,Ly-l,l)2j,Ly+2j-l = -— Y for 1 < i < [-£'j//2] (6.12) 
{TT,hc,Ly,Ly-l,l)2j,Ly+2j = {TT,hc,Ly,Ly-l,l)2j+l,Ly+2j = -— Y for 1 < J < [(1/2/ - l)/2] 

(6.13) 

(Tr,ftc,Lj„Lj,-l,l)Lj,+2j-l,2j-l = {TT,hc,Ly,Ly-l,l)Ly+2j-l,2j = ?/ " 1 for 1 < j < [ivj,/2] (6.14) 



{TT,hc,Ly,Ly-l,2)l,l — X 



(6.15) 
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{TT,hc,Ly,Ly-i,2)j,j ^i+x for 2<j<Ly-l 



(6.16) 



{TT,hc,Ly,Ly-l,2 



^Ly,Ly 



1 + X for Ly odd 
X for Lj, even 



(6.17) 



{TTM,Ly,Ly-i,2)2i,2i-,i = (7V,;,c,L„L,-i,i)2,+i,2,- = 1 for 1 < J < [(L, - l)/2] (6.18) 



(7T,;ic,Lj„Lj,-l,2)Lj,+2j-l,L^+2j-l — 1 for 1 < J < [i^jz/S] 
{TT,hc,Ly,Ly-l,2)Ly+2j,Ly+2j = y for 1 < J < [(Lj, - l)/2] 



(6.19) 



(6.20) 



{TT,hc,Ly,Ly-l,2)2j-l,Ly+2j-l - iTT,hc,Ly,Ly-l,2)2j,Ly+2j-l - — for 1 < J < [Ly/2] (6.21) 



{TT,hc,Ly,Ly-l,2)2j,Ly+2j = {TT,hc,Ly,Ly-l,2)2j+l,Ly+2j = ^ _ ^ for 1 < J < [(-Ly - l)/2] 

(6.22) 

{TT,hc,Ly,Ly-l,2)Ly+2j,2j = (Tr,ftc,Lj„Lj,-l,l)Lj,+2j,2j+l = 2/ " 1 for 1 < J < [(L^ - l)/2] 

(6.23) 

with all other elements equal to zero. 

For general x and y, TT,hc,Ly,Ly-i has rank equal to its dimension, 2Ly — 1. In the special 
case 7/ = (chromatic polynomial), the rank is reduced to (3Ly — l)/2 for odd Ly and 
to (3/2)Ly — 1 for even Ly (and may be reduced further for particular x). In the special 
case a: = (flow polynomial), the rank is reduced to Lj^ — 1 (and may be reduced further 
for particular y). These results are in accord with our previous general findings that [12] 
np{sq, Ly, d) = np(tri, Ly, d) and [15] np{sq, Ly, d) = np(hc, Ly, d). 

We illustrate these general formulas for the cases Ly = 2,3,4. For this purpose we use 
the abbreviations x'l = 1 + a;2 = 2 + a; and r^^ = y/(y — 1) as before. 



r,/ic,2,i,i 



xi 1 ra 

1 Xl Td I = lT,sq,2,l , 

V V y 



Tn 



T,hc,3,2,l 



/Xi 1 Td \ 

1 xi Vd l/v 

0x0 1/v 

V V y 

y 1 / 



T,hc,3,2,2 



X 1/v 
X 1/v 
1 

f X 1/f \ 

Xi 1 1/v Td 

1 Xi Vd 
1 
y I) y J 



(6.24) 



(6.25) 
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r,/ic,4,3,l 



( 


Xi 


1 













\ 




1 


Xi 








r-d 


1/v 













Xi 


1 





1/v 


rd 










1 


Xi 








Td 




V 


V 








y 

























1 





K 








V 


V 








y ) 



r,/ic,4,3,2 



fx 00 1/v 0\ 

xi 1 1/v rd 

1 r-d 1/v 

x 1/v 

1 

u u y 

yo 1 / 



(6.26) 



The Tx,hc,Ly,Ly-i, X — T,Z are obtained via eq. (2.8) from these auxihary matrices. For 
example, 



T,hc,2,l 



fl+x + x'^ 1+X {l+x)y{y-iy^ 
1 + x 1 + X + x'^ {1 + x)y{y — 1)~^ 

\ y-^ y 



(6.27) 



or equivalently, 



Z,hc,'^A 



v'^{3v'^ + 3vq + q'^) v^{2v + q) v'^{2v + q){l + v) 
v^{2v + q) v'^{3v'^ + 3vq + q'^) v\2v + q){l + v) 
v^ v^ v^{l+v) 



(6.28) 



We find that in general, for degree d = Ly — 1, one of the 2Ly — 1 eigenvalues, say that 
for j = 2Ly — 1, has a particularly simple expression: 

>^T,hC,Ly,d=Ly-l,j = 2Ly-l = (6.29) 

or equivalently, 

Xz,hc,Ly,d=Ly-i,j=2Ly-i = v'^''y-'\v + qf . (6.30) 
The expressions for the other 2Ly eigenvalues are, in general, more complicated. 
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VII. GENERAL RESULTS FOR CYCLIC SELF-DUAL SQUARE-LATTICE 

STRIPS 



In this section we consider the Potts model for famihes of self-dual strip graphs of the 
square lattice with fixed width Ly and arbitrarily great length L^., having periodic longitudi- 
nal boundary conditions, such that all vertices on one side of the strip, which we take to be 
the upper side are joined by edges to a single external vertex. A strip graph of this type will 
be denoted generically as Gd and, in more detail, as GoiLy x L-^). The family of Gd graphs 
is planar and self-dual. In general, the graph GoiLy x L^) has n= \V\ — L^Ly + 1 vertices, 
equal to the number of faces, /. One motivation for considering the Gd strip graphs is that 
they exhibit, for any Ly, the self-duality property of the infinite square lattice so that, by 
(1.13), 

T{GD,x,y)=TiGn,y,x) . (7.1) 

Equivalently, by (1.14), aside from a prefactor, the partition function is invariant under 
V q/v. This important duality property is used to derive the formula for the transition 
temperature of the ^-state Potts ferromagnet, which occurs at the self-dual point v — y^. 

In Ref. [17] we gave the general form for Z{GD,Ly x Lx,q,v) which, in our current 
notation with m given in terms of by eq. (1.4), is 

Ly + l 

Z{Gn, Ly X L„ q,v)^Yl '^^'^Tr[{Tz,G,,L,,dr] (7.2) 

d=l 

or equivalently 

Ly + l 

T{Gd, Ly X L„ X, y) = q^' ^ K^''^Tr[{TT,Gn,Ly,dT] (7-3) 

d=\ 

where the factor of q~^ arising from a prefactor of l/(a; — 1), and 



2d- 



J4) 



V 2 / 

'^d-l-3\d-, 



j=o \ J / 

The first few of these coefficients are 

«:(^)=?(g-2), k^'^ ^ q(q - l)(q - 3), k^^^ ^ q(q - 2)(q' - Aq + 2) . (7.5) 
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One has 

d 

«:W = J](g-s,,,) (7.6) 

k=l 

where 

s<i,fe = 4cos'(^) for A; = 1,2,. .d (7.7) 

and 

^{d) ^ ^(d) ^ ^{d-i) for (i= 1,2,... (7.8) 

The dimension dim{Tz,GD,Ly,d) — nziGo, Ly,d) was calculated in Ref. [17]; two relevant 
results are 

nziGD,Ly,Ly) ^2Ly , nziGD,Ly,Ly + l) ^1 . (7.9) 

Combining eq. (7.1), (7.3), and the fact that K^'^^ is a function of q which, in turn, is 

a symmetric function under x ^ y (cf. eq. (1.11)), it follows that the traces of the m'th 

powers of the transfer matrices themselves and their eigenvalues also have this symmetry: 

Tr[(TG,,L„.(x, y)r] = Tr[{TG,,^^,d{y, (7-10) 

^T,Go,Ly,di^^ y) = ^T,Go,Ly,diy^ ^) " (^-H) 

where we have explicitly indicated the functional dependence of the transfer matrices and 
eigenvalues on the variables x and y. 

A. Determinants 

We find 

detiTT,Go,Ly,d) = (xy)^-"-(«°'^--^'<^) (7.12) 
for the full range 1 < li < Lj^ + 1. 

B. Traces 

For the total trace, defined as 

Ly + l 

TT{TT,Gu,Ly) = E >^^'^TT{TT,Gu,Ly,d) (7.13) 

we find 

TT{TT,Go,Ly) = {xyf^ . (7.14) 
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C. Eigenvalue for d = Ly + 1 for A = Go 

For d = Ly + 1, the transfer matrix TT,GD,Ly,d reduces to a scalar, namely 



TT,GD,Ly,Ly+l - >^T,GD,Ly,Ly+l - 1 ■ (7-15) 

D. Transfer Matrix for d = Ly, A = Go 

The transfer matrix TT^GD,Ly,Ly has dimension 2Ly. We obtain the following general 
formula, which we write for the Tutte polynomial, since it has a somewhat simpler form. 

{TT,Gn,Ly,Ly)l,l^l + X (7.16) 

{TT,Go,Ly,Ly)j,j = 2 + x for > 2 and 2 < J < (7.17) 
iTT,Go,Ly,Lyh,j+i = iTT,G,,Ly,Lyh+i,j = 1 for l<J<Ly-l (7.18) 
{TT,Go,Ly,Lyhj ^ y for Ly + l<j<2Ly (7.19) 

{TT,Gn,Ly,Ly)j,Ly+J ^ for l<J<Ly (7.20) 

{TT,Gn,Ly,Ly)j+i,Ly+j = for 1<J<L,-1 (7.21) 

(Tt,g,,l„lJl,+,,, = 2/ - 1 for 1 < j < Ly (7.22) 
(Tt,g,,l„lJl,+,j+i = y - 1 for 1 < j < - 1 (7.23) 



with all other elements equal to zero. The elements of Tz^GD,Ly,Ly are given by the relation 
(2.27) so that Xz,Go,Ly,d = v^y\T,Go,Ly,d- 

For general x and y, TT^Go,Ly,Ly has rank equal to its dimension, 2Ly. In the special cases 
y — and x — 0, which yield the chromatic and flow polynomials, respectively, the rank is 
reduced to Ly and remains equal to Ly if both x and y are zero. 

We illustrate these general formulas for the cases Ly — 1,2,3: 
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Tt,Gd,1,1 - 



(7.24) 



V y 



/ Xi 1 Td \ 

1 X2 Td Td 

V V y 

V V y J 



(7.25) 



T,Gd,3,3 



I xi 1 \ 

1 X2 1 

I X2 Q Td Td 

V V y 

V V y 



(7.26) 



\ V y J 



Thus, in general, the upper left-hand submatrix has a main diagonal with the first entry equal 
to Xi and the other entries equal to X2. Adjacent to this main diagonal are two diagonals 
whose entries are f , and the rest of the submatrix is comprised of triangular regions filled 
with O's. In the upper right-hand submatrix the main diagonal and the adjacent diagonal 
below it have entries equal to r^, and the other entries arc zero. In the lower left-hand 
submatrix the main diagonal and the adjacent diagonal above it have entries equal to v, and 
the other entries are zero. And finally, in the right-hand lower submatrix the entries on the 
main diagonal are equal to y and the rest of this submatrix is made up of triangular regions 
of O's. For the lowest values Ly — 1,2, some of these parts, such as the triangular regions of 
zeros, are absent. 

Using eq. (2.27), it is straightforward to obtain the Tz,GD,Ly,Ly from these TT,Go,Ly,Ly 
matrices. For example. 



As corollaries of our general result for TT^GD,Ly,Ly we calculate the determinant and trace: 




(7.27) 



det{TT,Go,Ly,Ly) = {xy) 



Ly 



(7.28) 



which is the d — Ly special case of (7.12), and 



Tr{TT,Gn,Ly,Ly) = Ly{x + y)+2Ly-l. 



(7.29) 
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VIII. SOME ILLUSTRATIVE CALCULATIONS 



A. Square-Lattice Strip, Ly = 2 

The Potts model partition function Z{sq, Ly x m, BC, q, v) and Tutte polynomial 
T{sq, LyXm, BC, q, v) (BC = boundary conditions) were calculated for the cychc and Mobius 
strips of the square lattice of width L^^ = 2 in Ref. [11]. We express the results here in terms 
of transfer matrices Tx,sq,2,d, X — Z,T, via eqs. (1.6) and (1.17): 

( l + x + x^ xird\ TV 
Tt,s,,2,0^[^ ^ y j (8.1) 



_( q^ + ^qv + ^v' {q + 2v){l + v)\ 
^ZM-y ^3 v\l + v) ) ^^-^^ 

The matrices TT^sq,2,i and Tj'^sq,2,2 have been given above. The corresponding results for the 
Ly — 2 Mobius strip follow from our general formulas also given above. 

B. Triangular-Lattice Strip, Ly = 2 

We illustrate our results for the Ly = 2 cychc strip of the triangular lattice. We obtain 

^^■'"'^'°-l(y-i)(i + x + z/) y{i + y) ) ^^-^^ 

or equivalently, 

_ ( bv"^ + Aqv + q^ + v-^ {q + ?>v + v'^){l + v)\ , . 

lz,tri,2fi-y ^2(^^3^^^2) v\2 + v){l + v) ) ^^"^^ 

The matrices TT,tri,2,i and TT,tri,2,2 were given above. 

C. Honeycomb-lattice Strip, Ly = 2 

For the Ly = 2 cyclic strip of the honeycomb lattice we calculate 



r - T.U^' {l + x){l + x^)y{y-l) 

y- 1 y 



.5) 



or equivalently, 
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TzM,2,o - ^4(J + ^^) (8-6) 

hu = S-u^ + IQv^q + IQv^q^ + 5vq^ + (8.7) 

/ii2 = {2v + q){2v'' + 2vq + ?')(! + v) (8.8) 

The other matrices relevant for the strip were given above. 
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IX. APPENDIX 

In this appendix we include tables of the numbers Anz{sq, Ly,d) and nz{sq, Ly,d,±) 
discussed in the text. 

TABLE I. Table of Anz{sq, Ly, d) for strips of the square lattice. Blank entries are zero. The last entry 
for each value of Ly is the total number of partitions with self-reflection symmetry. 



Ly 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 




1 


1 


1 




















2 


2 


2 


1 


1 


















4 


3 


3 


3 


1 


1 
















8 


4 


6 


4 


4 


1 


1 














16 


5 


10 


10 


5 


5 


1 


1 












32 


6 


20 


15 


15 


6 


6 


1 


1 










64 


7 


35 


35 


21 


21 


7 


7 


1 


1 








128 


8 


70 


56 


56 


28 


28 


8 


8 


1 


1 






256 


9 


126 


126 


84 


84 


36 


36 


9 


9 


1 


1 




512 


10 


252 


210 


210 


120 


120 


45 


45 


10 


10 


1 


1 


1024 
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TABLE II. Tabic of numbers nz{sq, Ly,d,±) for strips of the square lattice. For each Ly value, the 
entries in the first and second lines are nz{sq, Ly, d, +) and nz{sq, Ly,d,—), respectively. Blank entries are 
zero. The last entry for each value of Ly is the total Nz,Ly,\- 



Ly {d, +) 


n 1 

u, + 


1,+ 




3, + 


A 1 

4,+ 


5,+ 


6,+ 


1 1 


Q 1 


n 1 
9,+ 


1 n 1 




(a,-) 


u, — 


1,- 


2, - 


3, - 


4, - 


5,- 


6, - 


7, - 


o 

8, — 


9, - 


10, - 


AT" 


2 


2 


2 


1 
























1 




















6 


3 


4 


6 


3 


1 




















1 


3 


2 


















20 


4 


10 


16 


12 


4 


1 


















4 


12 


8 


3 
















70 


5 


26 


50 


40 


20 


5 


1 
















16 


40 


35 


15 


4 














252 


6 


76 


156 


145 


80 


30 


6 


1 














56 


141 


130 


74 


24 


5 












924 


7 


232 


518 


511 


329 


140 


42 


7 


1 












197 


483 


490 


308 


133 


35 


6 










3432 


8 


750 


1744 


1848 


1288 


644 


224 


56 


8 


1 










680 


1688 


1792 


1260 


616 


216 


48 


7 








12870 


9 


2494 


6030 


6672 


5040 


2772 


1140 


336 


72 


9 


1 








2368 


5904 


6588 


4956 


2736 


1104 


327 


63 


8 






48620 


10 


8524 


21100 


24330 


19440 


11688 


5352 


1875 


480 


90 


10 


1 






8272 


20890 


24120 


19320 


11568 


5307 


1830 


470 


80 


9 




184756 
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